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This article has been published as a chapter in "Chemical Reactivity Theory: A 
Density Functional View", ed. P. K. Chattaraj (CRC Press, New York, 2009), ch. 
8, p. 105. In it, an overview of the relationship between time-dependent DFT and 
quantum hydrodynamics is presented, showing the role that Bohmian mechanics can 
play within the ab-initio methodology as both a numerical and an interpretative tool. 



I. INTRODUCTION 

Since the early days of Quantum Mechanics, the wave-function theory has proven to be 

very successful in describing many different quantum processes and phenomena. However, in 

many problems of Quantum Chemistry and Solid State Physics, where the dimensionality 

of the systems studied is relatively high, ab initio calculations of the structure of atoms, 

molecules, clusters and crystals, and their interactions are very often prohibitive. Hence, 

alternative formulations based on the direct use of the probability density, gathered under 

what is generally known as the density matrix theory [l| , were also developed since the very 

beginning of the new mechanics. The independent electron approximation or Thomas- Fermi 

model, and the Hartree and Hartree-Fock approaches are former statistical models developed 

in that direction 121 . These models can be considered direct predecessors of the more recent 

fl 

density functional theory (DFT) [3(, whose principles were established by Hohenberg, Kohn 
and Sham in the middle sixties. According to this theory, the fundamental physical 
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information about a many-body system is provided by single-particle densities in a three- 
dimensional space, which are obtained variationally within a time-independent framework. 
When compared with other previous formalisms, DFT presents two clear advantages: (i) it is 
able to treat many-body problems in a sufficiently accurate way and (ii) it is computationally 
simple. This explains why it is one of the most widely used theories to deal with electronic 
structure — the electronic ground-state energy as a function of the position of the atomic 
nuclei determines the structure of molecules and solids, providing at the same time the forces 
acting on the atomic nuclei when they are not at their equilibrium positions. At present, 
DFT is being used routinely to solve many problems in gas phase and condensed matter. 
Furthermore, it has made possible the development of accurate molecular dynamics schemes 
in which the forces are evaluated quantum-mechanically "on the fly". Nonetheless, DFT 
is a fundamental tool provided the systems studied are relatively large; for small systems 
standard methods based on the use of the wave function render quite accurate results Q]. 
Moreover, it is also worth stressing that all practical applications of DFT rely on essentially 



uncontrolled approximations [zj (e.g., the local density approximation 



the local spin- 



density approximation or generalized gradient approximations [8]), and therefore the validity 
of DFT is conditioned to its ability to provide results sufficiently close to the experimental 
data. 

As mentioned above, standard DFT is commonly applied to determine ground states in 
time-independent problems. Hence, reactive and non-reactive scatterings as well as atoms 
and molecules in laser fields have been out of the reach of the corresponding methodology. 
Nevertheless, though it is less known than the standard DFT, very interesting work in this 



direction can also be found in the literature 



y 



fig , where DFT is combined 



with quantum hydrodynamics [or quantum-fluid dynamics (QFD)] (QFD-DFT) in order to 
obtain a quantum theory of many-electron systems. In this case, the many-electron wave 
function is replaced by single-particle charge and current densities. The formal grounds of 
QFD-DFT rely on a set of hydrodynamical equations [lO, 11, 12|. It has the advantage of 
dealing with dynamical processes evolving in time in terms of single-particle time-dependent 
(TD) equations, as derived by different authors [14J]. Apart from QFD-DFT, there are other 



TD-DFT approaches based on similar grounds, such as the Floquet DFT 16|, [17| or the 



quantal DFT 



181 ] . Furthermore, we would like to note that TD-DFT does not necessarily 



require to pass through a QFD or QFD-like formulation in order to be applied 19] . As 
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happens with standard DFT, TD-DFT can also be started directly from the many-body 
TD Schrodinger equation, the density being then determined from solving a set of TD 
Schrodinger equations for single, non- interacting particles 121 ] . 

Although trajectories are not computed in QFD-DFT, it is clear that there is a strong 
connection between this approach and the trajectory or hydrodynamical picture of quantum 
mechanics j^, independently developed by Madelung 21], de Broglie 22] and Bohm 23], 
and which is also known as Bohmian mechanics. From the same hydrodynamical equations, 
information not only about the system configuration (DFT calculations) but also about its 
dynamics (quantum trajectories) is possible to obtain. This fact is better understood when 
the so-called quantum potential is considered, since it allows us to associate the probability 
density (calculated from DFT) with the quantum trajectories. Note that this potential is 
determined by the curvature of the probability density and, at the same time, it governs the 
behavior displayed by the quantum trajectories. Because of the interplay between probability 
density and quantum potential, the latter conveys fundamental physical information: it 
transmits the nonseparability contained in the probability density (or, equivalently, the wave 
function) to the particle dynamics. This prop erty , on the other hand, is connected with the 



inherent nonlocality of Quantum Mechanics 



241 ]. i. e., two distant parts of an entangled 



or nonfactorizable system will keep a strong correlation due to coherence exhibited by its 
quantum evolution. 

The purpose of this chapter is to show and discuss the connection between TD-DFT 
and Bohmian mechanics, as well as the sources of lack of accuracy in DFT, in general, 
regarding the problem of correlations within the Bohmian framework or, in other words, of 
entanglement. In order to be self-contained, a brief account of how DFT tackles the many- 
body problem with spin is given in Sec. [Til A short and simple introduction to TD-DFT 
and its quantum hydrodynamical version (QFD-DFT) is presented in Sec. IIHI The problem 
of the many-body wave function in Bohmian mechanics, as well as the fundamental grounds 
of this theory, are described and discussed in Sec. [TV] This chapter is conluded with a short 
final discussion in Sec. [V] 
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II. THE MANY-BODY PROBLEM IN STANDARD DFT 

There are many different physical and chemical systems of interest which are characterized 
by a relatively large number of degrees of freedom. However, in most of cases, the many- 
body problem can be reduced to calculations related to a sort of inhomogeneous gas, i.e., a 
set of interacting point-like particles which evolve quantum-mechanically under the action 
of a certain effective potential field. This is the typical DFT scenario, with an ensemble 
of N electrons in a nuclear or external potential representing the system of interest. DFT 
thus tries to provide an alternative approach to the exact, nonrelativistic iV-electron wave 
function \l/(riSi, . . . , t^sn), which satisfies the time-independent Schrodinger equation and 
where rjy and sn are the space and spin coordinates, respectively. Because the methodology 
based on DFT is easy and computationally efficient in its implementation, this theory is 
still enjoying an ever-increasing popularity within the Physics and Chemistry communities 
involved in many-body calculations. 

To understand the main idea behind DFT, consider the following. In the absence of 
magnetic fields, the many-electron Hamiltonian does not act on the electronic spin co- 
ordinates, and the antisymmetry and spin restrictions are directly imposed on the wave 
function \P(riSi, . . . , r^sx). Within the Born-Oppenheimer approximation, the energy of 
an iV-electron system with a fixed M-nuclei geometry R takes the following form in atomic 
units: 

E = -\\ , [V-V T 7 i(ri;r / 1 )]dr 1 + f v cxt (K, n) 7l (n) dr x + f l2 ^ V2 ^ dri dr 2 , (1) 

where 71(1*1) and 72(1*1, r 2 ) are the diagonal elements of 71 (r x ; t[) and 72(1*1, r 2 ; r' 1; r' 2 ), respec- 
tively, which represent the one-electron (or one-particle) density and the electron-electron (or 
two-particle) correlation function, commonly used in DFT and electronic structure theory. 
In principle, it might seem that all the information about the system necessary to evaluate 
the energy is contained in 71(1*1) and 72(1*1; r-2), and therefore one could forget about manip- 
ulating the wave function. However, in order to avoid unphysical results in the evaluation 
of the energy, it is still necessary to compute the wave function \&(riSi, . . . , t n sn) that gen- 
erates the correct 71(1*1) and 72(1*1; 1*2) densities. Equation (TjQ) is the starting point of DFT, 
which aims to replace both 7 i(ri;r' 1 ) and 72(1"!, 1:2) by p(r). If we are only interested in 
the system ground state, the Hohenberg-Kohn theorems state that the exact ground-state 
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total energy of any many-electron system is given by a universal, unknown functional of the 
one-electron density. However, only the second term of Eq. ([1]) is an explicit functional of 
p(r). The first term corresponds to the kinetic energy, which is a functional of the com- 
plete one-electron density function 7i(ri;r' 1 ). For A-electron systems the most important 
contribution to the electron-electron term comes from the classical electrostatic self-energy 
of the charge interaction, which is an explicit functional of the diagonal one-electron func- 
tion. The remaining contribution to the electron-electron term is still unknown. These 
two terms are a functional of the one-electron density, namely the "exchange-correlation" 
functional. Thus, it is possible to define a universal functional which is derivable from the 
one-electron density itself and with no reference to the external potential v ex t(R., r). Accord- 
ing to McWeeny 25| . we can reformulate the DFT by ensuring not only that a variational 
procedure leads to p(r) — which is derivable from a wave-function ^(r^si, . . . , rjySjv) (the 
so-called iV-representability problem) — , but also the wave function belongs to the totally 
irreducible representation of the spin permutation group A. From a mathematical point of 
view, the above proposition can be expressed (in atomic units) as 



E = min 

p— >7i derived from ty£A 



+ / tOR.rxMrOdrx + i f ll^l^hl^ ^ 



2 J r 12 

+ a • T/ 1 , TC , £corr[72(ri,r 2 )] >. (2) 

72 derived from W£A 

This equation shows the relationship between the one-electron function, 7 1 (r 1 ;r' 1 ), and the 
main part of the energy functional — the rest of the functional, which is the electron-electron 
repulsion, depends on 72(1*1, r2). The last term is also a functional of the one-electron 
density. In the new reformulation of DFT, the methodology is almost universally based 
on the Kohn-Sham approach and only differs in the particular way to model the unknown 
"exchange-correlation" term. 



III. TIME-DEPENDENT DENSITY FUNCTIONAL THEORY 



An extension of standard DFT is its TD version. This generalization is necessary when 
dealing with intrinsic TD phenomena. In addition, it preserves the appealing flavor of the 
classical approach to the theory of motion. 
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The rigorous foundation of the TD-DFT was started with the works by Bartolotti |1( 

n 

and Deb and Ghosh Jllfl . However, the proofs of the fundamental theorems were provided 
by Runge and Gross 12]. One of those theorems corresponds to a Hohenberg-Kohn-like 
theorem for the TD Schrodinger equation. The starting point for the derivation of the TD 
Kohn-Sham (KS) equations is the variational principle for the quantum mechanical action 
(along this Section, atomic units are also used): 



sm= /V(*)i 

J to 



\V(t))dt. (3) 



This variational principle is not based on the total energy because in TD systems the total 
energy is not conserved. The so-called Runge-Gross theorem then states that there exists 
a one-to-one mapping between the external potential (in general, TD), v ex t(r, t), and the 
electronic density, p(r,t), for many-body systems evolving from a fixed initial state, \l/(£o- 
Runge and Gross thus open the possibility of rigorously deriving the TD version of the Kohn- 
Sham equations. This procedure yields the TD Schrodinger equation for the Kohn-Sham 
electrons described by the orbitals <j>k(r,t), 

*^^ = # K sM)0 fc (r,t), (4) 

where the KS Hamiltonian is 

H KS {r,t) = -^V 2 + ^ KS [p(r,t)], (5) 

with a TD-KS effective potential, usually given by the sum of three terms, which account 
for external, classical electrostatic and exchange interactions. The latter is the source of all 
non-trivial, non-local, strongly correlated many-body effects. 

By construction, the exact TD density of the interacting system can then be calculated 
from a set of non-interacting, single-particle orbitals fulfilling the TD-KS equation (j4j) and 
reads 

N 

p(r,t) = ^|0 fc (r,t)| 2 . (6) 

fe=i 

Further analysis from the minimum action principle shows that the exchange (xc) potential 
is then the functional derivative of that quantity in terms of the density, 
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where S xc includes all non-trivial many-body parts of the action. The above equations 
provide the starting ground for further derivations of the theory. Thus, in addition to the 
TD-KS scheme, other variants have been proposed across the years, which include the TD 
spin-DFT, the TD current-DFT, the TD linear response DFT and the basis-set DFT j^. 
Each method has its range of applicability, but discussing them is out of the scope of this 
chapter. 

Here we focus on yet another implementation, the single-particle hydrodynamic approach 
or QFD-DFT, which it provides a natural link between DFT and Bohmian trajectories. The 
corresponding derivation is based on the realization that the density, p(r, t), and the current 
density, j(r,t) satisfy a coupled-set of "classical fluid", Navier-Stokes equations: 

dp(r,t) 



dt 
dj{r,t) 



-Vj(r,t), (8) 
PMM), (9) 



dt 
with 

P[p](r,t) = -i(*[p](m(T),H(t)]\*[p](t)), (10) 

being a functional of the density and with initial conditions p(r,t ) and j(r, t ). 

One can finally show that the above coupled equations translate into one single 
particle nonlinear differential equation for the hydro dynamical wave function $(r,t) = 
p(r, t) 1 / 2 e* ,s ( r '*) in terms of potential energy functionals: 



-^ + %H)*M) = ^, (ID 



with v e g[p] being given by 



r i ^-Eel-el . 5i?nu-el 5E XC 5T COTT SE ext 
VM = — + ^p- + ^p- + ^p- + ^p- (12) 

and where j(r, t) = p(r,t)v(r,t), with VS*(r,t) = v(r, t). For many particle systems, this is 
still an open problem (see next Section for a new discussion). In Eq. ( fl2l) . each term cor- 
responds, respectively, to the interelectronic repulsion energy, the Coulomb nuclear-electron 
attraction energy, the exchange and correlation energy, the nonclassical correction term to 
Weizsacker's kinetic energy, and the electron-external field interaction energy functionals. 
A judicious choice in the form of the above functionals yields surprisingly good results for 
selected applications. 
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As a simple mathematical approach to QFD-DFT, let us consider that the iV-electron 
system is described by the TD orbitals <j) k (r,t) when there is an external periodic, TD 
potential, for which we want to obtain the (TD) density p(r,t). These orbitals can be 
expressed in polar form, 

<f> k (r,t) = R k (r,t)e iS ^' t \ (13) 

where the amplitudes R k (r,t) and phases S k (r,t) are real functions of space and time, and 
the former are subject to the normalization condition 

i4(r,t)^(r,t)dr = 5 fc/ , (14) 

where f denotes the time-averaged integration over one period of time. The kinetic energy 
associated with this (noninteracting) iV-electron system reads now [h]] as 

i N r r 

T s [{R k ,S k }} t = --J2 / / {Rk(r,t)[V 2 R k (r,t)} - R k (r,t)[VS k (r,t)} 2 } dr. (15) 
fc=i J 

Similarly to the time-independent case, here we also assume the constraint that the sum of 
the squares of the R k gives the exact density p(r,t), i.e., 

N 

J>l(r,t)=p(r,t). (16) 
fc=i 

Moreover, we introduce an additional constraint: the conservation of the number of particles, 



k=l 



where j is the single-particle quantum density current vector. After minimizing Eq. (fl5l) with 
respect to the R k (which is subject to the previous constraints), we reach the Euler-Lagrange 
equation 

— -V 2 -Rfc + v eS R k = e k R k , (18) 

where v e s(r, t) and e k (r, t) are the Lagrange multiplier associated with the constraint defined 
in Eq. (Tl6|) and the conservation of the number of particles given by Eqs. (j!4p and fTTTj) . 
respectively. Moreover, e k (r,t) can be split up as a sum of two terms 

e k (r,t) = e k 0) +e k 1) (r,t). (19) 
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The quantity ejj, is a result of the normalization constraint, while ejj. are the Lagrange 
multipliers associated with the charge- current conservation defined by Eq. (ITT)) . On the 
other hand, if Eq. (fTSj) is divided by Rk we can reexpress the corresponding equation as 

Q k (r,t) + v cS (r,t) = e k (r,t) (20) 

where Q k is the so-called quantum potential associated with the state (ftk, 

Qk{r,t) = -l^. (21) 

Next, we minimize T s [{Rk, Sk}]t with respect to Sk to be subject to the constraint 

^ = -e k (r,t). (22) 

The resulting Euler-Lagrange equation is given by 

^ + V-(Rl VS k ) = 0. (23) 

The coupled equations, Eqs. (fT8l) and ff23|) . provide a means of determining the exact TD 
density of the system of interest. We note that, at the solution point, the current vector is 
given by 

TV 

j(r,t) = ^^(r,t)V4(r,t). (24) 
fc=i 

Note that, in the limit that the time- dependence is turned off, the TD-DFT approach 
correctly reduces to the usual time- independent DFT one, since VSfc vanishes, Eqs. (jTTj) . 
ff22|) and ff23l are identically satisfied, and Eq. fTl5|) will reduce to the time-independent 
kinetic energy of an iV-electron system. 



IV. BOHMIAN MECHANICS. A TRAJECTORY PICTURE OF QUANTUM 

MECHANICS 

A. Single-particle trajectories 



Apart from the operational, wave or action-based pictures of Quantum Mechanics pro- 
vided by Heisenberg, Schrodinger or Feynman, respectively, there is an additional, fully 
trajectory-based picture: Bohmian mechanics 20|, |23j. Within this picture, the standard 
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quantum formalism is understood in terms of trajectories defined by very specific motion 
rules. Although this formulation was independently formulated by Bohm, it gathers two 
former conceptual ideas: (1) the QFD picture proposed by Madelung, and (2) the pilot 
role assigned to the wave function, proposed by de Broglie. In this way, the time-evolution 
or dynamics of the system is described as an ideal quantum fluid with no viscosity; the 
evolution of this flow of identical particles is "guided" by the wave function. 

The Bohmian formalism follows straightforwardly from the Schrodinger one in the posi- 
tion representation after considering a change of variables, from the complex wave function 
field to the real fields (p, S) according to the transformation relation: 

tf(r,t) = R(r,t)e iS ^' t)/h , (25) 

with p = R 2 . Substituting this relation into the TD Schrodinger equation for a single particle 
of mass m, 



ih 



dV(T,t) 



2 



^ V2 + ^ 



tt(r,t), (26) 



dt 

and then separating the real and imaginary parts from the resulting expression, two real 
coupled equations are obtained: 



where 

cff ^ 2m R Am 



V z p 1 /V/; x 2 ' 



(28) 



P 2 V P 

is an effective potential resulting from the sum of the "classical" contribution, V, and the 
so-called quantum potential, Q, which depends on the quantum state via p — or, equivalently, 
on the instantaneous curvature of the wave function via R. Note that in the case V = v e s 
and \I> given as in the previous Section, one gets V e s = or V c s = Cfc(r, t) depending 
on whether we are considering the time-independent or the TD case, respectively. The 
action of the whole ensemble through the wave function on the particle motion can be seen 
as a dynamical manifestation of quantum nonlocality. Equation ( 127aj) is the continuity 
equation for the particle flow (or the probability density, from a conventional viewpoint) 
and (I27bl) is a generalized (quantum) Hamilton- Jacobi equation. As in classical mechanics, 
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the characteristics or solutions, S, of Eq. (I27bl) define the particle velocity field, 



(29) 



m 

from which the quantum trajectories are known. Uncertainty arises from the unpredictability 
in determining the particle initial conditions — distributed according to p(r, t = 0) 20J — , 
but not from the impossibility to know the actual (quantum) trajectory pursued during its 
evolution. 

An alternative way to obtain the quantum trajectories is by formulating Bohmian me- 
chanics as a Newtonian-like theory. Then, Eq. (129]) gives rise to a generalized Newton's 
second law, 

m^ = -W eff . (30) 

This formulation results very insightful; according to Eq. (!30|) . particles move under the 
action of an effective force, — W e ff, i.e., the nonlocal action of the quantum potential here 
is seen as the effect of a (nonlocal) quantum force. From a computational viewpoint, this 



formulation results very interesting in connection to quantum hydrodynamics [21|, [27| . Thus, 
Eqs. f l27|) can be reexpressed in terms of a continuity equation and a generalized Euler 
equation. As happens with classical fluids, here also two important concepts come into play: 



28 



which occur at nodal regions where the 



the quantum pressure and the quantum vortices 
velocity field is rotational. 

Since TD-DFT is being applied to scattering problems in its QFD version, two important 
consequences of the nonlocal nature of the quantum potential worth stressing in this regard. 
First, relevant quantum effects can be observed in regions where the classical interaction 
potential V becomes negligible and, more important, where p(r,t) ~ 0. This happens 
because quantum particles respond to the "shape" of but not to its "intensity", p(r,t) 
— notice that Q is scale-invariant under the multiplication of p(r,t) by a real constant. 
Second, quantum-mechanically the concept of asymptotic or free motion only holds locally. 
Following the classical definition for this motional regime, 

m^^O, (31) 
etc 

this means in Bohmian mechanics that VKfr ~ 0, i.e., the local curvature of the wave 
function has to be zero (apart from the classical-like requirement that V « 0). In scatter- 
ing experiments this condition is satisfied along the directions specified by the diffraction 
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channels 
forces. 



291 ]; in between, although V ~ 0, particles are still subject to strong quantum 



B. Bohmian trajectories describing many-body systems 

In the case of a many-body problem, the Bohmian mechanics for an iV-body dynamics 
follows straightforwardly from the one for a single system, but replacing Eq. ( |25i) by 

tt(r ls r 2 , . . . , r^; t) = i?(n, r 2 , . . . , r N ; t) e «(n,r a ,..,r Ar; t)/ft j (32) 

with p(ri, r 2 , . . . , r^; t) = R 2 (ri, r 2 , . . . , r^; t). If we are interested in the density of a single 
particle, we need to "trace" over the remaining N—l degrees of freedom in the corresponding 
density matrix (see Sec. IIVCI) . On the other hand, in order to know the specific trajectory 
pursued by the particle associated with the kth degree of freedom, we have to integrate the 
equation of motion 

v fc = , (33) 

m 

where Vfc = d/dvk- The velocity field is irrotational in nature except at nodal regions. 
Obviously, there will be as many equations of motion as degrees of freedom. Note that since 
each degree of freedom represents a particle that is interacting with the remaining N—l 
particles in the ensemble, the corresponding trajectory will be strongly influenced by the 
evolution of those other iV— 1 particles. This entanglement is patent through the quantum 
potential, which is given here as 

< 34 » 

fc=i 

where Q = Q(ri, r 2 , . . . , r^; t) is nonseparable and, therefore, strongly nonlocal. Note that 
this nonlocality arises from correlation among different degrees of freedom, which is different 
from the nonlocality that appears when considering symmetry properties of the wave func- 
tion, not described by the Schrodinger equation but by quantum statistics. In this sense, we 
can speak about two types of entanglement: symmetry and dynamics. The general iV-body 
wave function fl3"2"l) is entangled in both aspects. 

Now, if the many-body (electron) problem can be arranged in such a way that the many- 
body, nonseparable wave function is expressed in terms of a separable wave function which 
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depends on N single-particle wave functions (Hartree approximation), i.e., 

*(n, r 2 , . . . , r N ; t) = nf =1 ^ fc (r fe ; t) = U^ =1 R k (r k ; t) e iS ^' n , (35) 

then, in terms of trajectories, we find a set of uncoupled equations of motion, 

v fc = (36) 
m 

which will only depend implicitly (through v c s) on the other particles. Note that the fac- 
torization of the wave function implies that the quantum potential becomes a separable 
function of the N particle coordinates and time, 

fc2 N V72R N 

k=i ft k=i 

where each Q k measures the local curvature of the wave function associated with the ith 
orbital associated to the corresponding particle. Therefore, each degree of freedom can be 
studied separately from the rest (with the exception that we have to take into account 
the mean field created by the remaining N— 1 particles). Factorizability implies physical 
independence, statistical independence or, in other words, that particles obey Maxwell- 
Boltzmann statistics (they are distinguishable) and the associate wave function is, therefore, 
not entangled. 

In TD-DFT, the wave function is antisymmetrized and, therefore, nonfactorizable or 
entangled. However, as said above, it is not entangled from a dynamical point of view 
because the quantum forces originated from a nonseparable quantum potential like Eq. (Si 
are not taken into account . 



C. The reduced quantum trajectory approach 

In Sec. HVBt we have considered the problem of the reduced dynamics from a standard 
DFT approach, i.e., in terms of single-particle wave functions from which the (single-particle) 
probability density is obtained. However, one could also use an alternative description 
which arises from the field of decoherence. Here, in order to extract useful information 
about the system of interest, one usually computes its associated reduced density matrix by 
tracing the total density matrix, p t (the subscript t here indicates time-dependence), over the 
environment degrees of freedom. In the configuration representation and for an environment 
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constituted by N particles, the system reduced density matrix is obtained after integrating 
pi = \^>) tt (^>\ over the 3N environment degrees of freedom, {r^}^, 

p(r,r';t) = J (r, n, r 2 , . . . y N \m (t)) (t)\v\ n, r 2 , . . . r N ) d ri dr 2 • • • dr N . (38) 

The system (reduced) quantum density current can be derived from this expression, being 

j(r,£) = -Im[V r p(r,r';t)] , (39) 

771 r'=r 

which satisfies the continuity equation 

P + Vj = 0. (40) 

In Eq. ( I40p . p is the diagonal element (i.e., p = p(r, r;t)) of the reduced density matrix. 
Taking into account Eqs. (|39|) and (|40|) . now we define the velocity field, r, associated with 
the (reduced) system dynamics as 

J = pr, (41) 



which is analogous to the Bohmian velocity field. Now, from Eq. (j4ip . we define a new class 
of quantum trajectories as the solutions to the equation of motion 



. _ h Im[V r p(r, r'; t)\ 
m Re[p(r, r'; t)] 



(42) 



These new trajectories are the so-called reduced quantum trajectories 30], which are only 
explicitly related to the system reduced density matrix. The dynamics described by Eq. ( )42l) 
leads to the correct intensity (whose time-evolution is described by Eq. fj40|) ) when the 
statistics of a large number of particles is considered. Moreover, it is also straightforward to 
show that Eq. fj42|) reduces to the well-known expression for the velocity field in Bohmian 
mechanics when there is no interaction with the environment. 



V. FINAL DISCUSSION AND CONCLUSIONS 

Nowadays the success of DFT and TD-DFT is out of question in both the Physics and 
Chemistry communities. The numerical results obtained are most of cases in good agree- 
ment to those issued from experimental and other theoretical methods with a relative small 
computational effort. However, in this chapter, our goal has been to present the TD-DFT 
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from a Bohmian perspective and to analyze, from a conceptual level, some of the aspects 
which are deeply rooted in DFT. 

Working with a system of fermions, where the total wave function has to be antisym- 
metrized with respect to two-particle interchanges, it gives rise to the appearance of new 
quantum forces from the quantum potential which are not described by the DFT Hamilto- 
nian. The DFT wavefunction will be then nonfactorizable and, therefore, entangled from a 
symmetry point of view but not from a dynamical point of view. In this sense, as mentioned 
above, the effective potential V e g plays a fundamental role not only in the nonlocality of the 
theory, but in the so-called invertibility problem of the one-to-one mapping up to an additive 
time-dependent function between the density and v e g. In our opinion, the central theorems 
of TD-DFT should be written in terms of V e s instead of v e g, since the quantum potential 
is also state-dependent and a functional of the density. An infinite set of possible quantum 
potentials can be associated with the same physical situation and Schrodinger equation, and 
therefore the invertibility should be questioned. Moreover, for scattering problems, when 
v e s is negligible in the asymptotic region, the quantum potential can be still active and 
the time propagation should be extended much farther in order to obtain a good numerical 
convergence. 

In Bohmian mechanics, the way how the full problem is tackled in order to obtain opera- 
tional formulas can determine dramatically the final solution due to the context-dependence 
of this theory. More specifically, developing a Bohmian description within the many-body 
framework and then focusing on a particle is not equivalent to directly starting from the 
reduced density matrix or from the one-particle TD-DFT equation. Being well aware of the 
severe computational problems coming from the first and second approaches, we are still 
tempting to claim that those are the most natural ways to deal with a many-body problem 
in a Bohmian context. 
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